Spin and valley quantum Hall ferromagnetism in graphene 

A. F. Young 1 ", 1 C. R. Dean*, 2,3 L. Wang, 3 H. Ren, 1 P. Cadden-Zimansky, 1 
K. Watanabe, 4 T. Taniguchi, 4 J. Hone, 3 K. L. Shepard, 2 and P. Kim 1 

1 Department of Physics, Columbia University, New York, New York 10027, USA 
^ ^Department of Electrical Engineering, 

Columbia University, New York, New York 10027, USA 
3 Department of Mechanical Engineering, 
^ Columbia University, New York, New York 10027, USA 

* Advanced Materials Laboratory, National Institute for Materials Science, 
1-1 Namiki, Tsukuba, 305-0044, Japan 
^3 (Dated: January 20, 2012) 

t/3 

a 
a 

a 
o 
o 



> 



O 

(N 



1 



Electronic systems with multiple degenerate degrees of freedom can support 
a rich variety of broken symmetry states. In a graphene Landau level (LL), 
strong Coulomb interactions and the fourfold spin/valley degeneracy lead to an 
approximate SU(4) isospin symmetry. At partial filling, exchange interactions 
can break this symmetry, manifesting as additional Hall plateaus outside the 
normal integer sequence. Here we report the observation of a number of these 
quantum Hall isospin ferromagnetic (QHIFM) states, which we classify accord- 
ing to their real spin structure using tilted field magnetotrasport. The large 
activation gaps confirm the Coulomb origin of all the broken symmetry states, 
but the order depends strongly on LL index. In the high energy LLs, the Zeeman 
effect is the dominant aligning field, leading to real spin ferromagnets hosting 
Skyrmionic excitations at half filling, whereas in the 'relativistic' zero LL, lattice 
scale interactions drive the system to a density wave. 

The low energy effective theory of nearly neutral graphene describes two ffavors of mass- 
less Dirac quasiparticles centered on the two inequivalent corners of the Brillouin zone, 
termed valleys. In a high magnetic field, the valley degeneracy combines with the physical 
electron spin to produce four component Landau levels, leading to the anomalous graphene 
quantum Hall sequence 

4e 2 / 1\ 

where the LL index N is a nonnegative integer and the additional factor of 1/2 is related 
to the pseudospin winding number 1 . The spin/ valley degeneracy makes graphene a prime 
candidate for observing the rich physics of associated with multicomponent quantum Hall 
effects^. Graphene is exceptional as compared with its semiconductor counterparts due to 
a near-perfect energetic hierarchy (See Fig. [T]d). The energy scales characterizing cyclotron 
motion (En) and long range interparticle Coulomb interactions (Eq) — both of which reflect 
physics that is independent of spin or valley flavor — dwarf explicit spin and valley symme- 
try breaking effects. The combined four flavor degeneracy can therefore be thought of as 
that of a single SU(4) isospirP^. As in other multicomponent quantum Hall systems, ex- 
change interactions can drive the system through a ferromagnetic instability^, in which the 
order parameter corresponds to a finite polarization in a specific direction within the SU(4) 
isospin space. At integer fillings within a partially filled quartet LL, this order parameter 
is predicted to lead to a finite gap for charged excitations and a robust quantum Hall effect 
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for integers outside the sequence described in Eq. ([T]). The precise SU(4) polarization for 
given experimental conditions depends on the interplay between anisotropics arising from 
the Zeeman effect, lattice scale interactions, and disorder. All of these anisotropies are small 
and experimentally tunable, allowing for the possibility of a variety of distinct ground states 
across experimentally accessible ranges of filling factors, magnetic fields, and realizations of 
disorder. 

Previous studies have indeed reported the observation of QHE at several integer filling 
factors outside the normal sequence^^ however, the nature of the (presumably broken sym- 
metry) states leading to these plateaus remains a matter of intense theoretical debate^3HSl_ 
In the N=0 LL, most experimentaP^^H^ and theoretical^ work has focused on the 
strongly insulating behavior observed at z/=0, corresponding to half filling of the zero energy 
LL, which has no analog in conventional two dimensional electron systems. The insulating 
state has been described variously as a spin-polarized valley singlet, a valley polarized spin 
singlet, or a lattice scale spin density wave, but experimental resolution of this discrepancy 
has been hampered by the absence of any probe of the spin or valley order. Even less 
is known about the symmetry breaking at z/ == bl^ 15 1 17 1 18 1 23 ^ or throughout the iV 7^ 
LL d 8 | 9 | ll t ^H D ue to the anomalous structure of the N = LL, in which the valley quan- 
tum number corresponds to a real-space sublattice, the symmetry broken states for N = 
may not resemble those for N 7^ 0; however, limitations on sample quality and geometry in 
Si0 2 -supported and suspended devices, respectively, have precluded a comparative study. 

In this article, we address these issues by studying the thermal activation gaps, U A, asso- 
ciated with the broken symmetry IQHE states in graphene devices fabricated on hexagonal 
boron nitride (hBN) substrates^. This gap is associated with energy cost of the lowest 
lying charged excitations of the ground state. Owing to the atomic-scale confinement of the 
electronic wavefunctions to the plane of the graphene, all orbital effects related to electronic 
interactions depend only on the out of plane component of magnetic field (B±) while spins 
respond directly to the total magnetic field (Bt) independent of its direction. Tilted field 
measurements of the U A thus allow us to extract quantitative information about the net 
spin of charged excitations in the broken symmetry states. 

Fig. [T^i shows the evolution of the quantum Hall effect with magnetic field in a repre- 
sentative device (#165 thermal cycle 1). Symmetry breaking at v — ±1 is visible at fields 
of B± > 5 T, followed by the higher LLs at B± > 7 T; by 14 T (overlaid), p xx minima are 
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visible at all integer filling within the experimental range. In addition, an insulating state 
develops at v — (Fig. [TJ^) starting from 2-3 T, consistent with previous work on clean, 
suspended graphen d n * 37 l We find that all broken symmetry p xx minima, as well as the v = 
insulator, show simply activated temperature dependence over a wide range of fields (Fig. 

and d, insets; note that [T]l is data from a different device), allowing us to extract the 
energy gap U A as a function of perpendicular and total field. The exceptional quality of the 
devices studied here allows the observation of all integer filling broken symmetry states at 
magnetic fields of a few tesla, allowing the Zeeman energy to be tuned across a wide range 
in experimentally realizable magnetic fields. Using this technique, we explicitly demonstrate 
the dependence of the isospin ferromagnetic order, for fixed relative filling, on LL index N. 
For iV 7^ a dominant Zeeman anisotropy leads to spin polarized ground states at half filling 
and valley textured excitations at quarter filling, while in the N = LL, the situation is 
reversed: the v = insulator is shown to be unpolarized, and real spin textures form the 
most favorable charged excitations at from the fully polarized v — 1. 

We start our discussion with the half filled LLs. Figs. [2^, and d show activation gaps 
for the half-filled LLs, v = 0, -4, -8 and -12, as a function of B±. The perpendicular 
field dependence of U A for all N is qualitatively similar, following an approximately linear 
scaling with B±. We define the effective gyromagnetic ratio in perpendicular field from 
the slope of these cures, g± = /i^c^^A). We find that gu is enhanced with respect to 
the bare value g = 2. The Coulomb energy is the only scale in the system compatible 
with the measured energy gaps, which are much larger than might be expected from known 
single-particle effects. Moreover, g± decreases with increasing LL index N, consistent with 
exchange-driven quantum Hall ferromagnetisnP. 

Tilted field measurements reveal the uniqueness of the N = LL. For N ^ 0, the half filled 
gaps, ~ 4 A, -8 A, and _12 A, increase with Bf for fixed B± (see Fig. [2ji). This is consistent 
with the existence of real spin polarized states, in which excitations involve quasiparticles 
containing a net spin reversal relative to the ground state (and applied magnetic field). The 
activation gaps of such excitations consist of a direct Zeeman contribution from the reversal 
of spins against the external field, as well as an exchange contribution (Ax) arising from 
the spin reversal relative to adjacent (polarized) spins, 

A = A x (B ± )+g ol i B B T -T. (2) 
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In contrast, the v — resistance decreases (Fig |2p with increased Bt, an observation 
incompatible with the real spin polarized scenario 

for v = (J 22 I 25 I 28 I 31 I Instead; the data 

suggest that a spin unpolarized state, in which excitations contain a net spin aligned parallel 
(rather than antiparallel) to the applied field, underlies the insulating behavior at v = 0. 

Half filling of a fourfold degenerate graphene LL provides an ideal testing ground for the 
relative strength of the spin and valley anisotropies within the SU(4) isospin space. Because 
each cyclotron guiding center is doubly occupied, Pauli exclusion prevents the half-filled LL 
from fully polarizing in both spin and valley simultaneously. As a result, spin- and valley- 
polarizing tendencies necessarily compete, and the resulting ground state reflects the result 
of that competition. The fact that different order prevails at half filling for N = and N ^ 
even under identical experimental conditions (B± and Bt, which together fix the relative 
magnitude of the real spin anisotropy) suggests that the difference between LLs is intrinsic 
to graphene and originates in the valley sector. A likely origin lies with the unique structure 
of the ZLL wavefunctions: whereas for the N ^ LLs wavefunctions in a single valley are 
spread equally over the two real space sublattices, for the ZLL electrons in a single valley 
are localized on a single sublattice^. Long range interactions do not distinguish between 
such lattice scale orbital structural difference, but short range interactions do, potentially 
leading to different ground states in the N = and N LLs^Sl. At u = 0, the resulting 
interaction-induced valley anisotropies have been predicted to drive the system to one of a 
number of sublattice-ordered ground states I 20 I 21 I 23 I 24 I 27 I 29 I 3 ° I 33 E1] some of which are depicted 
in Fig. ^p. The experimental data presented here indicate that while the Zeeman effect 
wins the competition for the N ^ LLs, leading to spin polarized states at v — —12, —8, 
and —4, the valley anisotropies dominate the zero LL, leading to the formation of one of 
the possible lattice scale density waves portrayed in Fig. [2^ (ii), (iii), and (iv). The large 
size of the measure °A gap, and its insensitivity to in-plane fields, suggest that the valley 
anisotropies may be > 10 times stronger than their naive scale, {1/Ib) x Eq, in line with 
renormalization groupP^ and numerical calculations^. 

An additional notable feature of the experimental data is the linear scaling of °A on 
B±. The absence of a spin-polarized ground state at v — precludes linear-in-.Br Zeeman 
contributions to the excitation energy; meaning that the linear dependence must have an 
orbital origin. In contrast to 2D electronic systems with parabolic dispersion, in monolayer 
graphene both Eq and the cyclotron gaps scale as y/B±, a fact reflected in all available 
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theories of the v = insulator based on the low-energy Dirac model. The observed linear 
scaling thus points to the relevance of physics beyond the low energy theory. Motivating 
future work, quantitative data (See Supplementary Information) on the decrease of the gap 
in applied parallel field suggest that the spin ferromagnet, which is predicted 22 * 25 ^ to be an 
exact analog of the spin quantum Hall state^H may be experimentally accessible in the best 
samples at high tilt angles in accessible magnetic fields {Bt <45 T). 

Despite the role of the single-particle Zeeman effect in setting the order in the higher LLs, 
tilted field activation gaps demonstrate that the symmetry breaking can be thought of as 
being essentially spontaneous, with Zeeman functioning as a small aligning field. The gaps 
at half filling for N increase with total magnetic field (Fig. [2^) faster than might be 
expected for single spin flips, as reflected by the enhanced measured values of g\\ = fj,^ 8b t A 
(Fig. [2^, inset). In a Zeeman-dominated spin polarized state, charge transport occurs 
through the thermal activation of spin- reversed particle-hole pairs (Eq. [TJ. While exchange 
contributions to the energy gap (the first term in Eq. [I]) lead to g± > g^, this enhancement 
does not carry through to g» : changing Bt with B± fixed results in a measurement of the net 
spin of the excitation, and thus g\\ = go- In contrast, in exchange dominated spin polarized 
states it can be more energetically favorable to flip multiple spins smoothly in a Skyrmionic 
spin texture^ES, leading to a modified gap equation 

A = A X (B ± , K) + (2K + l)gofi B B T - T (3) 

where K > 0, the additional flipped spins per charged excitation, depends on the ratio 
Ez/Ec- The observed g» > 2 enhancements imply the presence of K > Skyrmionic charge 
carrying excitations at v — —4 and v = —8. 

The widely tunable field effect natural to graphene facilitates the study of transport 
phenomena over a wide range of density in the same sample, allowing a detailed analysis of 
the dependence of g\\ on LL index and B±. Our data is qualitatively in line with theoretical 
expectations for Skyrmion transport. The decreasing trend in g» with increasing LL index 
reflects the trend in exchange energy^. In contrast to two dimensional electronic systems 
with parabolic dispersion, where Skyrmions are expected only in the lowest energy LLp^HiS 
Skyrmions are theoretically possible for all |N|<3 in monolayer graphene^l for vanishing 
Ez- We observe signatures of Skyrmions up to N=2, consistent with the presence of a finite 
Zeeman effect. Measurements of g» at fixed filling for different values of B± in a single sample 
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(see Fi^3p) show the expected decrease in g\\ — for fixed Ec — with increasing E^: the Zeeman 
coupling punishes large Skyrmions, driving down K and, consequently, g». However, our 
data do not conform to the clean Skyrmion model, which predicts a universal scaling of the 
Skyrmion energy with the dimensionless parameter Ez/ Ec- The inset of figure |3^i shows the 
normalized in-plane dielectric constant, g\\/Ec, for three different values of B± as a function 
of Ez/Ec- While the two lower B± curves collapse onto one another, in accordance with 
theoretical expectations for the clean limit, the B±=25T data deviates strongly, implying 
the Skyrmions are larger than predicted 3 * 43 * 44 !. The failure of a universal scaling with Ez/Ec 
suggests that a second dimensionless parameter, for example related to the disorder strength 
(e.g., T/Ec) may be relevant for Skyrmion transport!^. 

We investigate the role of disorder in more detail by comparing several samples of varying 
quality. Overall, the values of g» observed in our graphene samples are lower than those 
measured — for similar values of Ez/Ec — in lower disorder GaAs quantum welLp^. Moreover, 
we find that g\\ is sample quality dependent. Figure^ shows tilted field dependence of ~ 4 A 
at B±=15T for several devices with different amounts of disorder. While higher disorder 
is correlated with smaller LL gaps, consistent with an increased LL broadening T, disorder 
also leads to a smaller measured g». In the most disordered samples measured here, we find 
that g\\ ~ 2, matching a previous measurement of SiCVsupported graphene^ and implying 
that in disordered samples, current is carried by single spin reversed particle-hole pairs. We 
speculate that Skyrmion size may by principally limited not only by the Zeeman effect but 
by the disorder landscape, which may tend to favor smaller Skyrmion size®'. The interplay 
between disorder and exchange may also contribute to the scaling of g\\ with B± (Fig. [3^,), 
which served to reduce the Skyrmion radius, r s k y oc (2K+ 1)£b 3 , in comparison to the length 
scale characterizing the disorder potential landscape. 

Like half filling, the phenomenology of a quarter filled graphene Landau level also shows 
markedly different behavior for N = and N ^ 0. At quarter filling the naive ground 
state is a fully polarized state in which a single spin-valley flavor is occupied. While spin 
is always polarized in the direction of the field, valley anisotropies are thought to lead to 
Ising or x-y type valley polarizations for N = and N ^ 0, respectivel y 15 * 17 ! Unlike at 
half-filling, spin and valley anisotropies do not compete with each other in the formation of 
the ground state at quarter filling: for a singly occupied cyclotron guiding center, there is 
no Pauli exclusion restriction on simultaneous spin- and valley polarization. The interplay 
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between spin- and valley- anisotropies does, however, contribute to the energetic of the 
excitation spectrum relevant for charge transport. In the case of a dominant Zeeman effect, 
for example, the low lying charged excitations are thought to consist of valley flip textures^ 
due to the high relative energetic cost of flipping real spins against the physical field. Most 
(although not all) activation measurements taken at quarter filling for N ^ (Figure 
and c) are consistent with this scenario, with the gaps independent of Bt for fixed B± to 
within experimental error. In contrast, gaps at v — — 1 increase with increasing Ez (Fig. 
|4p,), suggesting again that the Zeeman effect is not the dominant anisotropy in the zero 
energy LL. Tilted field dependence of the gaps at v — — 1 also show enhancement of g» over 
the bare value, consistent with the large strength of exchange in the zero LIP. 

While most data taken at half and quarter filling fit into the picture of QHIFM with 
LL index dependent anisotropy, we have observed several unexpected tilted field anoma- 
lies at both odd and even filling in the high LLs. For example, among samples for which 
~ 3 A was measured as a function of Ez, one (#129 (1)) shows reentrant behavior, with the 
gap collapsing with increasing Zeeman and then growing again as Ez is further increased. 
Qualitatively similar behavior involving gap collapse in parallel field was observed in other 
samples throughout the quarter filled N ^ LLs, (see Fig]^). Reentrant behavior was 
also observed at half filling in the higher LLs at low values of B± in one particularly high 
mobility sample (see Fig. [5)3). The origin of this unexpected behavior is not known, and 
the phenomenology is sample dependent: while the sample dependence makes the reen- 
trant behavior unlikely to be an intrinsic property of clean graphene, it tends to appear in 
cleaner samples, and disappears upon sample contamination. Among possible explanations 
are the unique susceptibility of the odd filling factors, to disorder^ES which can stabi- 
lize spatially inhomogenous ground states, or an effect of the hBN dielectric, which may 
induce a superlattice-like potential and sample dependent spatially inhomogeneous valley 
anisotropies. 

In summary, we have observed interaction driven spontaneously broken symmetry states 
at all integer filling factors in monolayer graphene. The picture that emerges is one of ex- 
change driven quantum Hall ferromagnetism within the combined spin- valley isopsin space^, 
where the leading anisotropies differ between the N=0 and N^O LL. Several questions re- 
main, however: first, the precise nature of the v — state remains elusive, and the linear Bj_ 
dependence is unaccounted for theoretically. Moreover, the absence of spin polarization at 
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v = and the presence of spin- reversed excitations at v = 1, need to be reconciled with the 
prevailing theoretical models of the fractional quantum Hall effect graphene N = LLl 49 * 50 !, 
all of which neglect the role of the valley anisotropics. Finally, little is understood about 
the structure and excitation spectrum of the odd filling states, and in particular the anoma- 
lous tilted field dependence observed in certain samples. We expect that the preliminary 
results presented here should motivate future work combining transport and surface science 
techniques, such as controlled absorption and scanned probe microscopy, to both elucidate 
the properties of these correlated states and, more generally, to use the graphene QHIFM 
as a model material platform for the systematic study of interacting systems in which all 
relevant experimental parameters, including disorder, can be tuned and probed in situ. 

I. METHODS 

We follow the fabrication method described in references IT2lfT3l to produce multiterminal 
Hall bar and van der Pauw graphene devices on hexagonal-Boron Nitride (h-BN) substrates. 
Most devices presented were etched into the Hall bar geometry using a short exposure to O2 
plasma; all samples were annealed in a 5% H 2 /95% Ar atmosphere to at 350° C for several 
hours to remove processing residues. Typical sample sized ranged from .5 x .5 to 5 x 5 /i m. 

Devices were measured in a sample-in- 4 He vapor variable temperature cryostat fitted 
with a mechanical sample rotation stage, mounted in the bore of a 35 T resistive magnet 
at the National High Magnetic Field Lab in Tallahassee, FL. Electrical measurements for 
v ^ were performed in the four point geometry using a 10-100 nA current bias. For the 
v = measurements were taken in the two terminal geometry using a 200 fiV excitation 
voltage. The numerous features present in a gate voltage trace at intermediate magnetic 
fields (B± < 25 T) allowed precise angle calibration, with B± determined to better than 
.5% accuracy. This was particularly important in the case of the v = state, where the 
dependence on B± is at least one order of magnitude stronger than that on Bt- 

R xx minima were determined by sweeping the gate voltage at fixed temperature. All 
p xx minima, as well as the resistance maximum at v — 0, obey an Arrhenius law, R xx ~ 
i?o ex p(^)- Gaps were determined by fitting to this formula over at least one decade 
of resistance when possible. Error bars are dominated by ambiguity in picking out the 
appropriate 'linear regime'. Plots for all gaps presented in the main text, including best fits, 
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are available in the Supplementary information (SI) 

Under normal conditions, graphene on h-BN samples are stable, and activation gaps can 
be measured repeatedly in the same sample over multiple cooldowns separated by weeks or 
months. However, rapidly warming the sample chamber causes outgassing from the cryostat 
walls, leading to the adsorption of debris on the graphene surface and higher disorder. This 
process, while difficult to control and not reliably reversible, allows us to exclude device- 
specific effects stemming, for example, from the interplay between the graphene electrons 
and the staggered lattice scale potential generated by the h-BN substrate. 
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FIG. 1: All integer quantum Hall effect in graphene on h-BN. a Landau fan from a 
monolayer graphene on h-BN device. Symmetry breaking of the LLs is visible from a few tesla, 
and at B± = 14T (superimposed) all integer filling factors feature minima in R xx . The color 
scale for —2 < v < 2 has been expanded by a factor of 7. b Energy scales. The cyclotron 
gaps (E N ~ ™pL±) a nd Coulomb energy (Ec 



jf- ) both characterize physics that does not 



distinguish between isospin flavors. The leading isospin anisotropies, which include the Zeeman 
effect {Ez = q^lbBt)-, lattice scale interactions (Esr ~ {a/ £b)Ec)"^ : , and disorder broadening(r), 
are at least one order of magnitude smaller. Here, vf — 10 8 cm/s is the Fermi velocity in graphene, 
a = 2.46A is the lattice constant, g$ = 2 is the bare gyromagnetic ratio, \ib is the Bohr magneton, 



and £b ~ 26 nm/yi?_|_ [Tesla] is the magnetic length. The disorder energy scale is extracted from 
the magnetic field dependence of the Shubnikov-de Haas oscillations at low magnetic field in one 
representative device (see Supplementary Information), c Development of the v = insulating 
state. Inset: temperature dependence showing the Arrhenius behavior of the insulating resistance, 
d Temperature dependence of the R xx minima in the symmetry broken IQHE regime. Inset: 
Arrhenius plots for i/ = 4asa function of magnetic field. 



15 




20 30 

B T 0) 



FIG. 2: a Activation gaps of half filled quartet Landau levels, a B± dependence of the v = 

gap, °A, for several devices. °A increases approximately linearly with applied B±, a feature not 

associated with any currently proposed theory for the v = 0. The dashed line indicates g± = 23. b 

Tilted field dependence of the resistance of the v = state. The resistance increases exponentially 

with field, consistent with a gapped state with °A oc B±. The resistance at fixed B± decreases for 

higher tilt angles, indicating a spin-unpolarized state, c Candidate QHIFM states for the v = 0. 

Our experiment rules out the spin ferromagnet, (i); all other states are marked by lattice scale 

spin (for the canted antiferromagnet, (ii)) or charge (for the charge density wave (iii) or Kekule 

distortion (iv)) order, d B± dependence of the half filled quartets for N ^ 0, v = —4, —8, —12. 

Like the v = 0, all gaps scale approximately linearly with B±, with enhanced g± factors that 

decrease with increasing LL index, e Unlike the v = state, all activation gaps measured for 

half-filled LLs with N ^ increase with Bt, indicating spin polarized states. For v = —4 and 

—8, the enhancement of gn indicates that charged excitations contain multiple flipped spins, f 

Schematic of charged excitations at half filling for N ^ 0. Excitations into the spin-reversed 

conduction band can take the form of single reversed spin particle hole pairs or smoothly varying 

Skyrmion-antiSkyrmion (S-aS) spin textures, depending on the strength of exchange interactions 

relative to disorder and the Zeeman energy. At B± = 15T in the samples studied in this work, the 

S-aS scenario prevails at u = —4 and —8, while charge at v = —12 is carried by single electron hole 
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pairs. 
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FIG. 3: Skyrmion transport at v = —4. a Perpendicular field dependence of gu at v = —4. 
The three curves show tilted field dependence of ~ 4 A for three different values of the Is- g\\, 
calculated from nearest neighbor and next-to-nearest neighbor finite differences of the curves in 
the main panel, is negatively correlated with Ez/Eq for fixed £b- For fixed Ez/Ec, g\\ is negatively 
correlated with £b, likely a combined effect of the increased exchange energy and decreased disorder 
parameter, T/Ec- Lines are guides for the eye. b Device dependence of _4 A for 2?j_=15T. Cleaner 
devices show both larger activation gaps (at all fillings) and larger for fixed is- This trend is 
observed both across devices and for the same device on different cooldowns, e.g. device #129. 
The inset shows the device dependence of gii , which decreases with increasing disorder (for which 
A(0=O)serves as a proxy variable). 
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FIG. 4: Transport phenomena at quarter filling, a v = — 1. Energy gaps increase with 
Bt, suggesting that excitations from the expected spin polarized, sublattice polarized stat c 1 15 * 1 ^^ 
involve real spin flips. Dependence on Bt appears to support enhancement of (solid lines are 
best linear fits to the data; dashed lines show gn = 2 for reference), (b-c) Tilt dependence of 
v = — 3 and v = — 5. Most samples show minimal dependence on Bt, consistent with theoretical 
predictions of valley-textured excitations 14 '. It is likely that in the absence of a single particle 
coupling field, these valley Skyrmions may be larg d 47 ! 4 ^ . Note reentrant behavior in sample #129 
(see Fig. pj). Dashed lines show g» = 2 for comparison. 
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FIG. 5: Reentrant QHE in tilted field in the higher LLs. (a) Gap collapse at odd filling. 
Gaps at quarter filling in the N ^ LLs collapse in applied in-plane field in selected samples. Gap 
collapse, or collapse and reemergence (cf. ( — 3) A in device #129 (1), Fig. [Z| is device dependent. 
In the two devices in which it was observed, the behavior did not survive a subsequent cooldown. 
b At very low perpendicular magnetic fields, an additional minimum is visible only at half filling. 
With the addition of an in-plane field, this minimum disappears, before reemerging again as Bt is 
increased further. Simultaneously, minima at quarter filling emerge with increasing Bt, indicating 
that they too can be spin active, increasing with Bt- 
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